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A summary of recent progress in Chiral Perturbation Theory (ChPT) 
at the two-loop level is given. A short introduction to ChPT is included, 
along with an explanation of the usefulness of developing ChPT for partially 
quenched QCD. Further, our recent work in partially quenched ChPT is 
reviewed, and a few comments are given on older work in mesonic ChPT 
at the two-loop level. In particular, we quote the present best values for 
the low-energy constants of the 0(p 4 ) chiral Lagrangian. 

1. Introduction 

This talk describes some of the work done at two-loop order in mesonic 
Chiral Perturbation Theory (ChPT) during the EURODAPHNE 1,11 and 
EURIDICE networks. A more extensive review of ChPT at this order can 
be found in Ref. pp. The aim of this talk is not to provide a full introduction 
or review of the two-loop work, but rather to concentrate on a few key issues. 
The outline of this talk is as follows: Sect. [2] gives a short introduction to 
ChPT, Sects. [3] and H] discuss partial quenching, why it is thought to be 
useful enough to warrant two-loop calculations in partially quenched ChPT 
(PQChPT), and why calculations in the PQChPT sector are much more 
challenging than those in standard ChPT. Our own work in PQChPT is 
also briefly reviewed, and finally Sect. [5] gives a brief summary of other 
existing work at the two-loop order in ChPT. 

2. Chiral Perturbation Theory 

Chiral Perturbation Theory was introduced in the papers by Weinberg, 
Gasser and Leutwyler O H] which build on earlier work within current 
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algebra and non-analytic higher order corrections. It should be noted that 
since a significant number of highly detailed introductions and lectures ex- 
ist [3(6], only a few of the main aspects are given here. The QCD Lagrangian 

£qcd = -\g^G^ + [*«z/9?z. + iQR^QR ~ m q (QrQl + QlQr)] (!) 

q=u,d,s 

is invariant under the chiral symmetry SU(3) L x SU (3)p when the masses of 
the up, down and strange quarks are set to zero. This symmetry is expected 
to be spontaneously broken by the quark-antiquark vacuum expectation 
value 

(??} = (QlQr + QrQl) ¥= (2) 

to the diagonal subgroup SU(3) V . Since this involves the spontaneous 
breaking of 8 generators of a global symmetry group, Goldstone's theorem 
requires the existence of 8 massless degrees of freedom and that their inter- 
actions vanish at zero momentum. ChPT is an effective field theory built 
on these eight massless particles which are identified with the pions, kaons 
and eta. This involves a long-distance expansion in momenta and quark 
masses. Such an expansion, called power counting, is possible because the 
interaction vanishes at zero momentum, and was worked out to all orders 
by Weinberg [2] . An example from ChPT is shown in Fig. [TJ 



(p 2 ) 2 (1/p 2 ) V = p 4 




Jd 4 p 




(p 2 ) (l/p 2 )p 4 = p 4 



Fig. 1. Illustration of power counting in ChPT. On the left are shown: The low- 
est order vertex, the meson propagator, a loop momentum integration and their 
respective powers of a generic momentum p. The examples on the right show two 
one-loop diagrams that count as 0(p 4 ) as compared to 0{p 2 ) for the lowest order. 

In order to perform the power counting expansion, higher order La- 
grangians need to be constructed. This has to be done in order to know 
the total number of parameters needed at a given order in the expansion. 
These parameters are referred to as low-energy constants (LECs). This clas- 
sification was done at 0(p 4 ) by Gasser and Leutwyler [31 H| and at 0(p 6 ) in 
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Ref. [7]. The number of LECs needed for the partially quenched case was 
determined in Ref. [H [9] . A summary of these results is given in Table [TJ 

Table 1. The number of LECs (physical + contact terms) at the various orders in 
mesonic ChPT. 





2 flavour 


3 flavour 


3+3 PQ 


o( P 2 ) 


F,B 2 


F ,B 2 


F ,B 2 


o( P A ) 


l\,h\ 7+3 


L\,H\ 10+2 


LI, HI 11+2 


0(p 6 ) 


c[ 53+4 


Cf 90+4 


K\ 112+3 



The main problem here is to determine a minimal set of LECs. No simple 
and straightforward procedure is known for the determination of such a set. 
Only if all of the LECs can be separately determined from "experiment" , 
or more generally from QCD Green's functions, can one be sure that the 
set of LECs is indeed minimal. Heat kernel methods allow to determine the 
divergence structure independently of Feynman diagram calculations. This, 
done at 0(p 4 ) in j3j[4j and at C(p 6 ) in [10], provides a very welcome check 
on actual two- loop calculations. 



3. Partially Quenched QCD 

One of the major applications of ChPT at present, and likely even more 
so in the future, is the extrapolation of lattice QCD results to the physical 
values of the light u, d quark masses. An overview of the many uses of 
ChPT in lattice QCD can be found in the recent lectures by Sharpe |11| . 
The main emphasis here is on the partially quenched aspect, which can be 
implemented in lattice gauge theory. In order to extract observables, one 
typically evaluates a correlator, e.g. a two-point correlator to obtain masses 
and decay constants. This correlator is evaluated in Euclidean space via the 
path integral (or functional integral) formalism: 



(0|(U 75 d)(x)(d7 5 «)(0)|0) = [dq][dq][dG] (u l5 d)(x)(d l5 u)(0) 



x exp 



(3) 



where the integral over the quarks and the anti-quarks can be performed 
and one obtains, schematically, 



[dq][dq][dG] (uj 5 d)(x)(dj 5 u)(0) exp 



i id y £ 



QCD 



OC 



(4) 
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[dG] exp 



valence 



QCD> 



~* sea 
gluonic 

where Ipc denotes the full Dirac operator with a specific gluon field config- 
uration but including the quark masses. The remaining integral over all the 
gluon degrees of freedom in Eq. @ is performed by importance sampling. 
The part labeled "valence" is connected to the external sources (hence the 
name), while the part labeled "sea" describes the effects of closed quark 
loops, not connected to any outside lines. Of course, gluons provide cou- 
plings between all these fermion lines if we look at the functional integral 
as a sum over Feynman diagrams. 

One major problem is that the determinant labeled "sea" is extremely 
CPU time consuming to evaluate. This has led to several approximations, 
the most drastic of which is the quenched approximation, whereby the sea 
contribution is completely neglected and only the gluonic and valence ones 
retained. "Unquenched" means in this respect that the sea determinant is 
included. However, the high CPU time requirements make it difficult to 
vary the quark masses very much in this part, and changing quark masses 
in the part labelled "valence" is indeed computationally much cheaper. In 
partially quenched simulations one thus varies the sea and valence quark 
masses independently of each other. There are good arguments in favour of 
this approach: 

• It is clearly superior to the quenched approximation. 

• More systematic studies of the input parameters may be performed. 

• It turns out that some quantities can be extracted from different ob- 
servables in this way. 

• Unlike the quenched approximation, it is continuously connected to 
the QCD case. 

However, a number of drawbacks need to be remembered: 

• It is not QCD as soon as quark masses are different in the valence and 
sea sectors. 

• It is not a bona fide Quantum Field Theory so the spin-statistics the- 
orem and unitarity relations are not satisfied. 

Especially the latter point might be important, since the derivation of ChPT 
from QCD relies heavily on unitarity. Nonetheless, one expects that at least 
close to the QCD case, PQQCD will have a low-energy effective theory 
similar to ChPT. 
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4. Partially Quenched ChPT at Two Loops 

The central problem in PQChPT is thus to mimic, within ChPT, the ef- 
fect of treating closed quark loops and quark lines differently. This problem 
is depicted schematically in Fig. [2j In some of the early calculations, the 
quark flow was inferred directly from the flavour flow in the ChPT vertices. 
An alternative approach, often referred to as the supersymmetric method, is 
more systematic. A series of bosonic ghost quarks with spin 1/2 is added to 
QCD. Due to the different statistics, these may cancel the effects of closed 
loops of valence quarks. This method is illustrated in Fig. [3l 

Mesons Quark Flow 

Valence 

XI. XI 

Fig. 2. The meson loop diagram on the left has different types of quark flow, both 
valence and sea quark as indicated on the right. 

The supersymmetric method was originally introduced for the quenched 
case in Refs. [T2j [131 E] an d later extended to the partially quenched case, 
see Refs. [151 EH E! an d references therein. Also, an instructive discussion 
about ChPT in the partially quenched sector is given in Ref. [16J. For 
practical purposes, the QCD chiral symmetry may be replaced by a graded 
symmetry which is (assumed to be) spontaneously broken to its diagonal 
subgroup: 

SU(n v + n s \n v ) L x SU(n v + n s \n v ) R SU(n v + n 8 \n v ) v , (5) 

where n v ,n s denote the number of valence and sea quark flavours. The 
"Goldstone bosons" now have both fermionic and bosonic character. A 
large amount of work exists at one- loop order, see the references in [11] . One 

Mesons Quark Flow Quark Flow Quark Flow Quark Flow 
Valence Valence Sea Ghost 

Fig. 3. The effect of adding ghost quarks to the different quark loops in the mesonic 
one loop diagram shown on the right. 
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stumbling block at two-loop order was the determination of the divergence 
structure and Lagrangians. Fortunately, it was realized [HJ [181 CES1 M that 
the work of |10] could be taken over by formally replacing traces by 
supertraces and the general number of flavours by the number of sea quarks. 

The final expressions at two-loop order are highly complex, which is due 
in part to the larger number of independent quark masses, but mainly to 
the peculiarities of the flavour neutral mesons. These do not have a simple 
pole structure but consist instead of several terms including a double pole: 



-iGfAk) 



1 



k 2 - 



■Xi 



■/?,. 



tlx 



+ 



R5 



+ 



D7T 

^"qii 



+ 



(k 2 -Xi) 2 k 2 -\i k 2 -x n k 2 



Xr, 



(6) 



where the various R coefficients consist of powers of ratios of differences 
of quark masses [131 US]- The presence of the many ratios leads to the 
extremely long expressions. The double pole in Eq. ([6]) is due to the fact 
that PQChPT is not a full field theory. Thus the valence quark loops cannot 
be resummed to all orders, as shown in Fig. HI 



Full 




Quenched 



O CO 



Fig. 4. Schematic resummation of quark loops. Imagine gluons everywhere on top 
of the quark lines drawn. Each bubble gives a lowest order meson propagator. The 
full resummation exponentiates and leads to a single pole. For valence quarks only 
the loops shown at the bottom appear, leading to double poles in the propagator. 



A large amount of work exists in PQChPT at two loops. The first calcu- 
lation was the meson mass for the case of three sea quarks, with a common 
valence mass and a common sea quark mass [Sj. Since then the decay con- 
stants [T5] and masses [S] have been fully worked out. These quantities are 
also fully known for two sea quark flavours [TP]. The formulas are rather 
lengthy and the number of parameters is also quite high. Analytical pro- 
grams have therefore been posted by the authors on the website [6]. Methods 
of dealing with the large number of parameters have been discussed exten- 
sively in Refs. [19] and [9] for the two- and three-flavour cases, respectively. 
It should be emphasized again that the LECs of unquenched ChPT, and 
thus the full QCD results, are related to the partially quenched LECs in 
a simple way via the Cayley-Hamilton relations of Ref. [7]- More recent 
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work has focused on the neutral mass sector of PQChPT. It was shown in 
Ref. [16J how the residue T> of the double pole 

— iZV 

G U k ) = {k 2 _ M 2J2 + ( ? ) 

can be measured on the lattice and used at C(p 4 ) to extract Ly, which is 
relevant for the r] mass. In Ref. [20j the self-energy resummation was worked 
out to all orders, and it was shown explicitly how to obtain the double pole 
and the structure of the full propagator from the one-particle irreducible 
diagrams. It was also shown that all 0(p 6 ) parameters relevant for the r\ 
mass can be extracted using this method. The most recent work has focused 
on the inclusion of dynamical photons in the partially quenched theory [21] . 



5. Standard ChPT at Two Loops 

The existing work at two-loop order in mesonic ChPT is very briefly 
reviewed here. A much more extensive review may be found in Ref. pQ. 
The oldest two-loop work in ChPT made use of dispersive techniques to 
calculate the non-analytical dependence on the kinematical variables. This 
was done numerically [22] and analytically [23] for the pion vector and 
scalar form factors, and fully analytically for irir scattering in Ref. |24| . 
The first full two-loop calculations appeared somewhat later in the two- 
flavour sector, with 77 — > 7r°7r° |25] and 77 — > 7r + 7r~, F w and [26J. The 
process 77 — > -kit was recently recalculated in Ref. [27] . With -kit scattering 
[28], pion vector and scalar form factors [29] and the radiative decay of the 
pion [30J, most processes of interest have now been worked out. 

The earliest three-flavour work, on the vector two-point functions, was 
by Golowich and Kambor [31 j . extended to all flavour cases in Refs. [32 1 1331 
134] . The first calculations with proper two-loop integrals were of the meson 
masses and decay constants, in Refs. [351 E3] and [36], including isospin 
violation. All scalar two point functions [371 EH] and vacuum expectation 
values |39] are also known. More recent work covers the electromagnetic 
form factors [40 [ |4"I] , K i3 |4CH I42j , and scalar form- factors [13] . Processes 
with more external legs include K ell4 [39] , kaon radiative decay [H] , tttt [UJ 
and nK [JB] scattering. The first results at finite volume have appeared 
recently [471135] . 

A major problem in phenomenological applications of ChPT at two loops 
is to find enough experimental inputs to determine the 0(p & ) parameters. 
In practice most of them have to be estimated, which is typically done along 
the lines of Ref. [19] by saturating the LECs by resonance exchange. While 
this can be done at various levels of sophistication, most phenomenolog- 
ical applications have used a fairly simple extension of Ref. [49] . see e.g. 

Refs. [251 m\ EE El m S3 M ■ 
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Table 2. The fitted optimal values of the L\ at C(p 4 ) and 0(p 6 ), and convergence 
behviour of several quantities for the different fits [36l [46] . 





fit 10, 0(p 6 ) 


fit 10, 0(p 4 ) 


fit D 




i n 3 tx 




u.oo 


n aa 




W 3 L r 2 


0.73 ±0.12 


1.59 


0.69 




W 3 L r 3 


-2.53 ±0.37 


-2.91 


-2.33 




10 3 L r 4 


= 


= 


= 0.2 




10 3 L r 5 


0.97 ±0.11 


1.46 


0.88 




W 3 L r 6 


= 


= 


= 




10 3 L r 7 


-0.31 ±0.14 


-0.49 


-0.28 




10 3 L r 8 


0.60 ±0.18 


1.00 


0.54 




2B m/ml 


0.736 


0.991 


0.958 


ml 


0(jp% 0(p 6 ) 


0.006, 0.258 


0.009, = 


-0.091, 0.133 


m\ 


: 0(p 4 ), 0(p 6 ) 


0.007, 0.306 


0.075, = 


-0.096, 0.201 


m\ 


C?(p 4 ), 0(p 6 ) 


-0.052, 0.318 


0.013, = 


-0.151, 0.197 


m u /m d 


0.45 ± 0.05 


0.52 


0.50 


F K 

Frr 


F [MeV] 


87.7 


81.1 


80.4 


0(p 4 ), 0(p 6 ) 


0.169, 0.051 


0.22, = 


0.159, 0.061 



Most phenomenological applications rely on the work of Refs. [391 136j . 
The fitting method and the inputs used are described in detail in Ref. [36], 
and the results can be found in Table El in the columns labeled "fit 10" . 
These used the (at that time) most recent data of the BNL E865 experiment 
as the main K i4 input. The change compared to a fit at 0(p 4 ) is also given 
in Table [2j These fits assume that the 1/N C suppressed LECs L\ and Lg 
vanish at the scale \i = 0.77 GeV. On the other hand, "Fit D" of Ref. 06] 
uses all the same inputs as "fit 10" , in addition to the dispersive results on 
7T7T and ttK scattering from Refs. [50J and [51J. The convergence properties 
of some quantities are also given in Table El However, an update of the 
fit is in order, with the new experimental results on K i4 and an improved 
treatment of the 0(p 6 ) constants along the lines of Ref. [52"1 153"] . 

6. Conclusions 

ChPT at two-loop order is by now a very well developed field, where a 
large number of two- and three-flavour calculations have been performed. 
The use of the partially quenched results will hopefully allow for many of 
the 0(p e ) LECs to be determined from Lattice QCD, thus removing a major 
stumbling block in phenomenological applications. 
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